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Initial-switching refers to the way in which the decay of an initially confined state begins, as the 
barrier isolating it from the exterior is relaxed. We study these effects in the context of Longhi's 
version of the Fano- Anderson model. Most authors assume the sudden approximation where the 
coupling is turned on instantaneously. We consider a finite rise time T, both numerically and 
analytically. When the coupling is ramped up linearly over a switching time T, we show that the 
asymptotic survival amplitude acquires a phase T and is modulated by a factor (sin T)/T. Several 
other results relating to the solution of the model are obtained. All site amplitudes have the same 
decay constant during the exponential decay regime. In the asymptotic regime, the amplitude and 
phase of decay oscillations depend on the initial-switching profile, but the period does not. 
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I. INTRODUCTION 

Quantal effects in the time evolution of decaying sys- 
tems have been studied since the beginning of quantum 
mechanics. Gamow's explanation of alpha-decay was 
among the first successful applications of quantum the- 
ory to radioactive nuclei The theory showed that 
decay was exponential in the accessible range of times, 
and related the lifetime to the nuclear charge Z and de- 
cay energy. Khalfin [2j was the first to realize that, when 
the energy spectrum of the system is bounded from be- 
low, exponential decay cannot persist as t — > oo. Even- 
tually it has to be replaced by an asymptotic regime 
with a slower rate of decrease. Explicit solution of the 
Schrodinger equation for simple models of a particle es- 
caping from a confining potential verified that in the 
asymptotic regime, the survival (or non-escape) proba- 
bility shows power-law decrease Q . Under rather general 
assumptions, this power-law can be shown to be cx t~ 3 for 
finite range potentials [J]. Some additional aspects are 
discussed in [a, [6[, and references therein. Exponential 
decay is observed in a wide variety of physical systems. 
The predicted slower decay rate at long times has proven 
difficult to verify, but it is a universal feature present in 
theoretical models. This applies not only to a particle 
trapped by a potential barrier, but also to the decay of 
a discrete state coupled to a continuum @, § - 

It is less well-known that when the confining bar- 
rier, or the coupling between discrete state and contin- 
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uum, is time-dependent in the initial stages of decay, so- 
called "memory effects" affect the post-exponential sur- 
vival probability of the system. These effects were pro- 
posed as a means of enhancing the survival probability, to 
make it feasible to observe the post-exponential regime 
experimentally. In particular, Robinson Q considered 
adiabatic switching-on of the system over an infinite time 
period, to clarify the role of switching time versus obser- 
vation time. 

Tight binding models have been widely used in con- 
densed matter for many years [3, [H| to describe tightly 
bound bands in crystals. More recently they have been 
applied to multilayered semiconductor heterostructures 
, to the dynamics of Bloch oscillations [HI and to the 
study of properties of Bose condensates in a lattice [l4| • 
They are particularly well suited when the unit cell of the 
periodic struture leads to a well confined lowest energy 
state. 

In this paper, we examine memory effects for finite 
switching times. We use a variant of the Fano- Anderson 
model [la, lla LLZl proposed by Longhi [l8| to study de- 
viations from exponential decay at short times, and the 
Zeno effect [l9j]. His model can be viewed as the decay 
of a discrete state coupled to a single-band tight-binding 
continuum. But, more interestingly from an experimen- 
tal point of view, he presented it as a one dimensional 
semi-infinite periodic chain of discrete sites. The cou- 
pling between adjacent sites, g(\n)(n + 1| + \n + l)(n\), 
is described by a universal hopping constant, g, and the 
site energies are taken all to be equal. For two sites iso- 
lated from the rest, izU/(2g) would be the half-period, 
the time for a complete transfer of occupation between 
the pair. Relative isolation of the first site is obtained 
by fixing its coupling parameter g± = gA to be different 
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from the rest. Under these conditions, when at t = only 
site 1 is occupied, Longhi obtained the survival probabil- 
ity, its lifetime for exponential decay and the asymptotic 
power- law [H, H3] ■ The connection to models of a con- 
fined particle, like alpha-decay, views site 1 as the well, 
the ratio A = gi/g as the effect of a confining barrier, 
and the semi- infinite chain of sites in position space as 
a single-band continuum. Fig. [TJ illustrates the system 
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FIG. 1: Semi-infinite array of identical ID wells coupled to a 
"defect" well separated by a thicker barrier. 

we have in mind. The barrier separating site 1 from the 
rest of a semi-infinite array is thicker than the others, to 
reduce tunnelling, which in a tight binding model implies 
g± < g. A thick enough barrier effectively makes g\ = 0, 
removing site 1 from consideration. 

An added interest of the model is that it has a classi- 
cal electro-magnetic analogue, in the form of a waveguide 
array [2l|, [22j. Fig. [2] describes one such semi-infinte 
array with an initial "defect". Fig. 1(a) of j22j shows 
an experimental realization of such an array. Consider 
an ideal semi-infinite array of long parallel waveguides. 
The variation with z, (the longitudinal distance along 
the guide) of the electric field in the n'th guide, replaces 
the variation with time of the site amplitudes. Both sys- 
tems obey a set of coupled mode equations formally the 
same as those for the time variation of the site ampli- 
tudes in the Fano- Anderson model. Recent experiments, 
using scanning tunnelling microscopy with sensitive near- 
field probes, have allowed Longhi and collaborators [23| 
to measure the evanescent fields along the waveguides, 
and provide a convincing quantitative demonstration of 
the classical-quantal analogy for the case of a periodic 
parallel array. In addition the classical analogue of the 
quantum Zeno effect was verified 24] . The validity of the 
Fano- Anderson model to describe such effects has been 
quantitatively shown by these authors [H, [24| . 

In the present work, we study the problem of initial 
switching, by letting g%(t)/g = A v (t) depend on time, 
during a finite initial time period < t < T, following 
which A is constant. In Longhi's optical analogue, this 




FIG. 2: A segment of a semi-infinite array of tunnel-coupled 
optical waveguides, made by etching. The z-axis is parallel to 
the guides, and zt at the TOP of the drawing. A logarithmic 
dependence on z of the separation leads to linear dependence 
of the coupling. 



could be implemented by having many parallel waveg- 
uides coupled to a first curved waveguide, with its separa- 
tion distance decreasing smoothly with z from an initially 
large value (where tunnelling is practically negligible), to 
a constant value beyond z = Zt- The final separation dis- 
tance should be chosen to give well-expressed exponential 
decay. Longhi et al. (26j have demonstrated the experi- 
mental feasibility of a similar device (see their Fig. 1). 
Our work differs from that of Robinson , who consid- 
ered a smooth switching function increasing adiabatically 
over — oo < t < oo, so his A#(t) saturates only asymptot- 
ically. Finite switching times have also been considered 
in the context of threshold effects in pulsed laser photion- 
ization; see Chapter 6 of [25| for a review. In Section II 
we discuss solutions of Longhi's model in the exponential 
and asymptotic regimes. In Section III we derive ana- 
lytic expressions for the time evolution during the period 
of varying coupling strength. Specific calculations of the 
effects are reported in Section IV. Verification of these 
predictions should be possible with present experimental 
capabilities. Verification of the predicted memory effects 
at long times may require further development of exper- 
imental techniques. 



II. DECAY FOR ARBITRARY INITIAL 
CONDITIONS 

To begin, we study the time evolution of a semi-infinite 
tight binding system that initially has a single site n = q 
occupied at t = 0. For convenience, we will use dimen- 
sionless units. Energies will be measured in units of g 
(the hopping parameter of the periodic lattice), and time 
in units of h/g. Wth the site energies set to zero, the 



3 



Hamiltonian is 



H = -A„(t) 



|1><2| + |2)(1| 



E 



\n)(n + l\ + \n+l)(n\ 



(1) 



We assume a positive coupling parameter < A < 1. 
Longhi [HI, H(| derived analytic expressions for the oc- 
cupation amplitude c n (t) , focussed primarily on the case 
that only the first site q = 1 is occupied at time zero, and 
constant A. Fig. Oshows that the choice A = 0.3 gives a 
survival probability P(t) = |ci(i)| 2 , which decays expo- 
nentially for a considerable range of times, and afterward 
the post-exponential decay consists of oscillations whose 
envelope follows the t~ 3 law. This will be our example 
of reference in the remainder of the article. 
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FIG. 3: Survival probability P(t), for constant A = 0.3, in the 
region of changeover from exponential to asymptotic decay. 
Here and in the subsequent figures, time t ( dimensionless) is 
in units %/g. 

In this paper we suppose that during an initial period 
< t < T, the coupling parameter A„ increases contin- 
uously from to A and afterwards remains constant. In 
Section III we will study a specific choice for the initial 
time-dependence, that we call linear rise 

A v (t) = A e (t) = t A/T , t<T 

= A , t > T. (2) 

For this model we will derive explicit expressions for the 
survival amplitude and other quantities, valid to lowest 
order in A. Exponential decay over a sizeable time scale 
requires that A be small compared to 1, so linear terms 
should suffice to establish the importance of such effects. 

Initially we focus on general properties that apply to 
any scheme for the initial switching-on of A. Given any 
profile A v (t), the site occupation amplitudes, c n (t), n — 
1, 2, • • • can be determined by numerical solution of the 
time-dependent Schrodinger equation (TDSE) , for which 
we use the method described in [23. Alternatively, in 
Sec. Ill we use the time evolution operator to derive 



approximate analytic results. In either approach, the end 
result is a set of occupation amplitudes at t — T which in 
this Section we assume to be known, and that we write 
in vector form 



0&l $2 



(3) 



with Y^Li \^q\ 2 = 1 an d where by T we denote the trans- 
pose. Since we will always assume that after time T, 
A v (t > T) = A remains constant, it is convenient to 
reset the zero of time by writing t = t — T . In the re- 
mainder of this section we study the time evolution of 
a system with these evolved initial conditions, at t = 0, 
given by Eq. ([3]), and constant A. The further evolution 
is governed by the TDSE 
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with 



|$t)) = (c 1 (t),c 2 (t),..,c n (t),..) T . 

The linearity in time of the TDSE implies that 

00 

q=l 



(4) 



(5) 



(6) 



where the c n , q (t) are the site amplitudes of a state, 
\ip q (t)}, which at t = has unit real amplitude c 9) g(Q) = 1 
on the q-th site, and vanishes on all others. 

|^=o)) T = (o,o,- ■• ,i,rv-0 

\Mi)) T = (ci, q (t),c 2 , q (t),---c n , q (t),---) (7) 

Stationary solutions. To determine these c„ j9 (t), we 
first find stationary solutions of the coupled equations, 
and then propagate linear combinations of them in time. 
Writing \ip(t)) = exp(— iEt)\u), the components of \u) 
satisfy 

= 2 cos (j> ui — A U2 

= — Aui + 2cos(/>U2 — 1*3 

= — u n -\ + 2 cos 4>u n — u n +x , n>3 (8) 

where we have written E = —2 cos0. The solutio is most 
easily obtained by expressing u n in terms of incoming 
and outgoing waves with amplitudes denoted X and O: 



ui = A{X-0)/{2 cos0) 
Ami = 2 cos <j) 1*2 — 1/3. 



,71 > 2 



(9) 



Inserting the first and second lines into the third, gives 



O 2 cos 6 e 



A 2 



2 cos - A 2 



Y2i<t> 



r(<f>) , (10) 



where a 2 = 1 — A 2 is the decoupling parameter. Since the 
normalization of the states \u) is arbitrary at this stage, 
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we may define the amplitudes O and X as the numerator 
and denominator of Eq. (fTU|) . Then the u n in Eq. © 
take the compact form 



u n((f ) ) — 2i ( sin n0 + a 2 sin(n — 2) 
ui((j)) — 2iAsin0 , 



n > 2 



(11) 



(Notice that if line one is specialized to n = 1, we obtain 
A times iti (</>). Therefore one can avoid treating ri = 1 
as a special case, if we simply remember to remove this 
extra factor.) 

Time dependent solutions: Following 23] we write 
the c n (t) as weighted integrals over the stationary so- 
lutions, with a weight function F{4>) (still to be deter- 
mined) expressed as a Fourier series 



F{cj>) 



^ oo 

r — 7- > Tjv e~ 
2ttA ^ 



<•„(/! = / # e- iEi F(<j>) 

7T 



(12) 



The parameters T/v are determined by the initial condi- 
tions, as we now explain. In what follows we will often 
use properties of Bessel functions which can be found in 
Chapter 11 of Arfken [28|; to avoid tiresome repetition, 
we will cite it no further. For convenience we introduce 
symbols J q (2t) = i q J q (2t), which have the symmetry 
J q = J- q , and satisfy the recurrence relation 



\Jn—\ \Jn-\-\ — ..\Jn- 
it 



(13) 



This allows us to write 



Ip.N — — 



sin»m e 



i(2t cos c/>-N4>) 



= t N -PJ N _ p (2t)-i N +PJ N+p (2t) 

= J N _ p (2t)-J N+p (2t) (14) 

= JN~p — JN-p+2 + JN-p+2 — JN~p+i + 

• • • + 3n+p-i — JN+p 
(N + p + l-2m) 



E 



it 



-JN+p+l~2m ■ (15) 



In the limit of large argument 2t ^> N + p, J p jv ap- 
proaches its asymptotic form 



r (°0 

1 p.N 



E 



{N+p- 



2m) 



N+p-2m 



V^/J 3 / 2 



cos [2/5- (N +p + 1 - 2m+ l/2)?r/2] 



^ (N+p+l-2m) 



i n +p x 



cos[2i -{N + p + l + 1/2)tt/2] 

N+p-2 P 

2 E ( N + p + 1 ~ 2m ) x 



rn—1 

cos[2i- (N + p- 1/2)tt/2] 



,-AT+p-2 



TV p COS 



f 3 / 2 

Using Eqs. dHJ», OH) and (H 

= 4 E Tj 



2f 



N+p- 



1\ 7T 



2/ 2 



(16) 



'- JV 



JV=1 
.2 



{JN-n{2t) — jN+n(2t)) 



d(t) 



+ a 2 (jw-n+2(2<) - J" W +ri-2(2f)) 

n > 2 , while 
51 T JV (Jjv-i(2*) - Jw+i(2t)) 



JV=1 



EiV 
T N -~J N (2t) 

N=l lt 



E Tn 

N=l 



N - 1 jJ N {2t\ 



(17) 



so that at t = 0: 



Ac„(0) = (T„ + a 2 T„_2) 
ci(0) = Tx. 



71 > 2 



(18) 



Note that Tq — 0, since only the sites with positive in- 
dices exist. Given initial conditions for the c n (t), as e.g. 
in Eq. ([5]), one determines the Tjv, N — 1, 2, .. from the 
recurrence relation Eq. (jTHJ) . This completes construc- 
tion of the time-dependent solution. 

Up to now, our results are valid for all times t > 0. 
In the next two subsections we develop approximations 
valid in the asymptotic and exponential regimes which 
are illustrated in Fig. 1. 



A. Asymptotic approximation 

Since Jn{21)/1 cx (t)~ 3 / 2 when t — > oo, one sees that 
for long times |ci(i)| 2 cx t~ 3 for any set of initial occupa- 
tions, as stated in the Introduction. Using Eqs. (Til))) and 
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(fTTf we find 



Cl (t)^c{ a \t) 

v^L NT "' N -' 



(~ 7r 

cos 2t - N 

V 2 4 



^(2Z + l)r 2Z +i J sin (2t- j 



sZ = 



E 2 ^ 



\Z = 



2/ | COS \2t — 



(19) 



where for later convenience we have isolated the contri- 
butions from even and odd order Tjv. Providing that the 
Tjy are real, (which they will be,) this is equivalent to 
separating the real and imaginary parts. 

For the asymptotic behaviour of the c„, n > 2, we use 
Eqs. IHl), JIB]) and (T7]) to arrive at 



ci a) (t) 

E 



■n-l 



< 3 / 2 A 
N ~ 1 N T N cos ( 2t 



n + a 2 (n- 2) x 



N=l 



l\ 7T 

^+"-2 2 



■ra-l 



n + a 2 (n — 2) 



^(2Z + 1)T 2;+1 J cos (2* - (n + 1/2) |) 
z=o / 

i (E^ 2Z T2 ' ) cos ( 2i ~- («- 1/2) | 



. (20) 



for n > 2. Providing that 2i >> rt7r, we can replace the 
cosines in the last lines by J n (2t) and J n -i(2t) respec- 
tively, (making appropriate adjustments to the prefac- 
tor). 

Now we particularise to cases where the initial ampli- 
tude is confined to a single site labelled q. By combining 
these results, we can finally write the general solution for 
arbitrary initial occupations. 

Case 1: Initial amplitude only on site 1. At t = : 
Ci(0) = 1, c„(0) = 0, n > 2. It is then immediate from 
Eq. dHJ) that T 2[ = for all I and that T 2 i+i = (-) l a 21 . 
Inserting these in Eq. (flU]) one finds 



i 



^a 2Z 



z=o 

1 - a 2 . / ~ 7T\ . , 

— 7 ^sin (2f (21) 

0rf 3 / 2 (1 + a 2 ) 2 V 4/ v y 



in agreement with Eq. 8 of 18] . Hereafter we will call this 
case the "sudden approximation" because it corresponds 
to the abrupt switching-on of the coupling from A„ = 
to A„ = A at t = T = 0. 



Case 2: Initial amplitude on another odd site. Con- 
sider now c„(0) = 5 n>q , q > 1 odd. Then T 2 i = 0, I = 
0, 1, .. and 

T\ — T 3 = ■ ■ ■ = T q -2 = , T q = A 
T q+ 2i = (-) l a 2l A , (22) 

so that Eq. (fH)]) gives 



A q(l + a 2 ) - 2a 2 
T^W (1 + a 2 ) 2 



(2t-J) (23) 



Note that for q = 1 this result agrees with Eq. (|2ip. if 
we remove the explicit A, in agreement with our earlier 
remark. The meaning of the factor A is this: When site 
1 is initially empty, the wave must cross the barrier at 
least once before its decay may begin. 

Case 3: Initial amplitude on an even site. When 
c„(0) = 5 n>g , q even, we find T 2 i+i = 0, I = 0, 1, • • • 
and 

T = T 2 = ■ ■ ■ = T q _ 2 = , T q = A 
T q+2l = (-) l a 2l A. (24) 

From Eq. ((TSJ) 

(a) (f] -»A g(l + a 2 )-2a 2 
Cl ' q[) V^~ 3/2 (1 + a 2 ) 2 4 



cos 2*- 7 (25) 



General initial conditions. We now combine the pre- 
vious results to obtain the most general solution. Insert- 
ing Eqs. (USD and (US]) into © we find 



c^\t) 



1 



^(l + a 2 ) 2 f 3 / 2 
So sin (2t - t) - iSe cos (it - 



(26) 



with the amplitudes summed over odd (even) indices 

oo 

S o(e) = ^[9(1 + " 2 ) -2a 2 ], (27) 

q— odd(even) 

where we have defined $1 = $i and $ q = ^ g A, q — 
2,3,---. 

In principle, the phases of the amplitudes "d q at time 
zero are arbitrary, so that So and S e are neither purely 
real or imaginary. Then c^°' of Eq. (|26|) is a complex 
number and will have a node only accidentally. But in 
the case we wish to discuss, it is assumed that at some 
time in the past the system was prepared with only site 
1 occupied and c\(t = 0) = 1. During the switching on 
process, the system evolved to build up the amplitudes 
in sites 2, 3, • • • . In the special case of sudden switching, 
T = 0, we can apply Eq. (f20|) with t — t to see that the 
resulting phase of c q (t) is the quantities Tn being 

real. Using a perturbation expansion, we will see that the 
same phases are generated by the time-evolution operator 
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in Section III (see discussion following Eq. (|6T]) ). As a 
result, the amplitudes on odd- numbered sites are real and 
those on even numbered sites, pure imaginary. Hence S a 
and X e = —iS e are both real, allowing us to write 



» IT 



it) 

A 
4> 



v/SF(l + a 2 ) 2 t 3 /2 
arctan(X e /S ). 



„4sin {2t- j + 



(28) 



This simplification ensures that the survival probability 
will still have nodes, and the period of the asymptotic 
oscillations is not changed by the initial switching pro- 
cess. Compared to the sudden approximation Eq. (|2ip . 
there is an extra phase $, and the asymptotic amplitude 
involves A in place of A. To have $ = would require 
that the even-sites are empty at t = 0. 



B. Exponential decay regime 

Again, we first consider the situation where at t = 0, a 
single site q is occupied, and use these as building blocks 
for the general result to be given at the end. Inserting 
the expressions Eqs. ((211) or flU} into Eq. ([T2"]) 



1 



2tiz1 z 2 



(29) 



where we have set z = exp(i^) and the subindcx q refers 



to the initial condition, c„(0) = 
has poles at the origin and z p 
(fl"2l) becomes: 



$5,n- We note that F q (<j)) 
±ia. In terms of z, Eq. 



Cfc(i) 



1 



u k (z) 



zi- 



z 2 + a A 



(30) 



the path of integration being the unit circle. As shown 
by Longhi [3, Eqs. (9) to (12), to extract the 
exponentially-decaying contribution to Ck {t) one may use 
the residue theorem to evaluate the integral, retaining 
only the contribution from the Gamow pole at z g = —ia. 
The conjugate pole corresponds to exponential growth, 
which is not the physical process, and the poles at the 
origin give compensating contributions, leaving only a 
small remainder, at least in the survival probability. The 
Gamow pole contributes 



>) (t) = e -7*72 



(31) 



where the superscript (e) means exponential decay, and 
1 - a 



1 
2 



A2 

a 



(32) 



This is a remarkable result: it shows that all the ini- 
tial conditions of the form c„(0) = S„ tq lead to the same 



lifetime, r = 1/7, for the exponentially decaying con- 
tribution to |cfc, 9 (i)| 2 , no matter which site q is initially 
occupied. 

We develop the above result by evaluating Uk(z g ). 
From Eq. dTTJ) , for k > 2, 



Uk(Zg) 



•k 1 - " 4 

1 — ■ 

or 



and therefore: 



k+q A 2 (l + a 2 ) n % 
k^ b > - 1 2 a k +i 



(33) 



k > 2. (34) 



When k = 1, the same calculation has a factor 1/A in 
comparison: 



u\(z g ) = -i A 



1 + a 2 



and 



^ = _. g+1 AO + o^) 



2a?^ 



(35) 



(36) 



The sudden approximation, Case 1, with ci(0) = 1, re- 
quires a separate discussion since T21 = and T21+1 = 
(— ) 2 a 2 ', without the factor A of Eq. (|22]) . One now finds 



2ttA(z 2 + a 2 



and 



2a k+1 



when k > 2, while 



-7*72 



(37) 



(38) 



(39) 



This last result agrees with Eq. (11) of fl8j ]. 

Finally, we write the amplitudes for general initial oc- 
cupations. Inserting the above results in Eq. ^ gives 




-A(1 + « 2 )(X:^] e^ 2 , (40) 



when k > 2. Clearly the exponential decay contribution 
has a delay T, manifest in the appearance of t = t — T in 
these expressions. Furthermore, one sees that the decay 
constant 7 depends neither on the initial conditions nor 
on the site considered. 



III. INITIAL SWITCHING 

We now address our main objective, which is time evo- 
lution when the coupling between site 1 and site 2, A v (t), 



varies continuously with time. To do so we introduce a 
simple diagrammatic technique which describes visually 
the action of the time evolution operator. To introduce 
these diagrams we first apply them to the oft-studied case 
of "sudden switching": A v (t > 0) = A =constant, which 
will serve as a reference for results with variable A. 
The time evolution operator can be written as 



U(t) = U m{t) with 



m=0 



Um{t) = (-*)' 



dt r , 



dt m -i ■ ■ 



I 2 



x H(t m )H{t m _ 1 )---H(t 1 ) 



(41) 



and in the special case of a time independent H, this 
reduces to U(t) = exp(-iHt) = Em=oH*- ff )' T 7( m! )- 
Suppose we start with the state |1 > with only site 1 
occupied initially. Then 



U(t)\l) 



' Cl (t)> 

|*(t)) = I C2(*) 



(42) 




FIG. 4: Diagrams which begin at site 1, with m = 1 and 
2 lines. The meaning of the various symbols is explained in 
text. 



A. Diagrammatic method and sudden switching 

We now develop a diagrammatic method to calculate 
the time evolution operator, first in the context of sudden 
switching where the exact solution is already available. 
It will be extended to linear switching in the next subsec- 
tion. To begin we compute U m (t)\l) for increasing values 
of m, assuming A is constant. For m — 0, Uo(t) — 1, and 
therefore 



When 7ti = 1 

E/i(t)|l) 



U (t)\l) = (43) 



-% I dti H |1) = -itH |1) 
lo 

it A 1 2). (44) 



The first graph in Fig. Q] describes this process: the cir- 
cles represent the sites, the time evolution is downwards 
and the lines describe hopping between adjacent sites. 
The m — 2 contribution is 



u 2 (t) = H) 2 
u 2 (t)\i. 



dt x H 2 = 



2! 



o Jo 
(A 2 |l) + A|3)) 



2! 



t 2 H 2 



(45) 



and the two corresponding diagrams are in the lower row 
of Fig. 2J Continuation to higher m is straightforward. 
The rules to generate the contributions to U m (t)\l) are: 
1) each line joining 1 and 2 contributes a factor —A, 
whereas the other lines contribute a factor —1 each; 2) 
if the diagram has m lines: multiply by (— it) m /(m\). 




FIG. 5: Three line diagrams which begin at site 1. 



Combining these two, the net result is to multiply by 
(ii) m /(m\). The final state is the last vertex at bottom 
of the diagram. As a further example consider the m = 3 
contributions: applying the rules to the three graphs of 
Fig. [5l one can immediately write 



tMt)|i> 



3! 



A 3 |2)+A|2) + A|4) 



(46) 



We now focus on contributions to c\ (t) , to all orders in m. 
The relevant diagrams will be those having |1) as both 
initial and final state. Since this implies an even number 
of lines, we write m = 21 with I = 1,2, • • • . The simplest 
of these diagrams, I = 1, is already included in Fig. 0] 
The two diagrams with / = 2 are shown in Fig. [6l In the 
special case A = 1, each diagram with the same number 
of lines gives an identical contribution to c\. Applying 
the general rules, the desired series is 



d(*,A = l) 



oo 

E 

1=0 



(-) 



j-21 

(20! 21 



(47) 
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Eq. HI]) tells us that when A = 1 (a 2 = 0) and Case 1 
initial conditions, 





FIG. 6: Four line diagrams which begin and end at site 1. 

where N^f* 1 is the number of distinct diagrams with 21 
lines having |1) as both initial and final state. To de- 
termine it, we rewrite Eq. (fTT|) for the Case 1 initial 
conditions and A = 1: it simplifies to 



Cl (t,A = l) = 



Ji(2t) 



t 



E 

1=0 



I 4-21 



ll(l + l)V 



(48) 



using the Taylor series for Ji(2t). Comparing Eqs. (j4"T|) 
and |H8"1). it follows that 



(20! 



Z!(/ + l)! 



(49) 



Now that the iVj 1- * 1 are known, we can return to the more 
general case A ^ 1 and consider those diagrams that give 
contributions of order A 2 to c\ (t) . These begin and end 
at |1), but never have |1) as an intermediate state. The 
second diagram in Fig. [5] is an example. It is now easy 
to write down their contribution: 



,(2) 



I 



1 (20 



^ 2 E- 

i 



(-)' t 



l +21 



-(Ai-2)(i\y 



.(50) 



Here we have noted that the number of diagrams with 
2(1 — 1) lines and beginning and ending at |2) without 
|1) as intermediate states, equals Afc^n- Since we are 
interested in values of A < 1, we will not write down 
expressions for contributions of order A 4 or higher. This 
supposes that processes with more than two tunneling 
events between sites 1 and 2 can be neglected during the 
short initial-switching time interval. 

To find the amplitude for site 2 we begin by considering 
again the special case A = 1. The simplest diagrams that 
connect |1) to |2) are shown in Fig. 0](in the upper row) 
and Fig. [5] (the first two). They have an odd number of 
lines, m = 21 + 1, I = 0, 1, 2, • • • . Again all the diagrams 
with the same number of lines give equal contributions. 
Thus 



c 2 (i,A = l)=£ 



(it) 



21+1 



1=0 



(21 + 1)1 



(51) 



'•->"• A 1 ' = ¥ J 2(2<) = 2z<£ AJlL (52) 



t 



and therefore 



_ 2(2/ + !)! 
21+1 ~ 11(1 + 2)1 



(53) 



We now return to the case A 1, and determine the 
contribution to c 2 linear in A. The simplest diagrams 
are those in the top row of Fig. 0] and in Fig. More 
generally: we include all diagrams that begin with a link 
from 1 to 2 and are completed with lines that go from 
2 to 2 without having 1 as an intermediate state. All 
such diagrams have an odd number n = 21 + 1 of lines, 
I = 0, 1, • • • and their contributions are 



(*) = A E 



21+1 



1=0 



(21 + 1)} 



Ni 



*aE(-)' 



21+1 



1=0 



(21 + 1) 11 (/ + !)!' 



(54) 



where we have used that the number of the above sub- 
diagrams, beginning and ending at 2, equals N^ 1 ■ 
Finally, to determine the amplitudes for site n we proceed 
similarly: For constant A = 1, 



.(^=i) = E(-)'#t 



n-l+21 



1=0 



(n-l + 21)! "~ 1+2( 



(55) 



and comparison to Eq. (|17[) with Case 1 initial conditions 
c„(t,A = l) = z"- 1 ~J n (2t) 

j-n-l+21 

- in ~ ln J2(-y^rTnl (56) 



1=0 



l\ (n + l)\ 



shows that 



JV n-l+2( 



i(n- 1 + 21)1 
l\(n + l)l ' 



(57) 



which generalizes Eqs. (|49|) and (|53|) . Finally, for con- 
stant A ^ 1 we obtain the contribution to c„ linear in 
A 



f— ' n — 1 



n - 1 (it) 7 



1=0 



n-l + 2ll\(n-l + l)\ 



(58) 



B. Linear switching 

When A(t) depends on time, the time evolution opera- 
tor is given by Eq. (|4ip and the term of order m involves 
m consecutive integrations over intermediate times. In 
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computing the contribution linear in A to c n (i), the time- 
dependent term will appear only in the first line of the 
diagram. Thus for m = 1 



S e to lowest order in A: The derivation is entirely 
similar to that of S , and leads to 



CM* 2 )|1) = i I 2 dtx A v (h) |2) = f(t 2 )\2) (59) 
Jo 

and we construct the contributions from the diagrams 
corresponding to higher orders by integrating f(t) the 
required number of times. 



dtr 



dt ,- 



dt 2 f(t 2 ) (60) 



This l m replaces the t m /ml found when A is constant. 
For the "linear switching" case, Eq. ((2]) and t < T, the 
integral above gives (since f(t) = t 2 /2), l m = t' m+1 / (m+ 
1)!. The modified rule for diagrams corresponding to 
linear switching is to replace t m /m\ by t m+1 /(m + 1)! 
and also add the factor 1/T from Eq. ([2]). Thus, instead 
of Eq. (|58p . taking all relevant diagrams into account we 
find 



it) 



T 5^ 



j n ~ l+2s ^.n-\-2s 



(n-l) 



T (n + 2s)(n - 1 + 2s)s\ {n - 1 



(61) 



for t < T. Note that the phase of this amplitude is 
i n . This remains true if we include terms of higher 
order in A. As discussed following Eq. (|45|) . to connect 
site 1 with site n, the corresponding diagrams must have 
n — 1 + 2s lines with s = 1,2, .... Therefore the phase 
will be i™- 1 + 2;5 and terms with odd (even) n will be real 
(imaginary) . 

Our next task is to apply these amplitudes to evaluate 
the sums defined in Eqs. ()27[1 . to obtain approximations 
to the asymptotic amplitudes: 

S to lowest order in A: Since 

(1 + a 2 )q - 2c? = 2(q - 1) + A 2 {q - 2) ~ 2(q - 1), (62) 



SP = 2AJ2 (2;-l)41 } (T) 



.A2 

' T 



5»<(2;-l) 2 5> 

1=1 s=0 



(I + s)(2(l + s) - l)s!(2Z + s - 1)! 
.A 2 



i — sin T, 
T 



(65) 



where the sum is again detailed in the Appendix. To pic- 
ture how the functions S and S e behave, we computed 
them by numerically solving the TDSE to determine the 
q = 1, 2, — , and then computed S a and S e (also nu- 



merically) using Eqs. (|27| . Fig. [7] shows the curves 
thus obtained for three values of A. As expected, for 
the smallest A the curves are very close to the predic- 
tions of Eqs. (EU) and flBSJ), but even for A = 0.9, the 
approximation is qualitatively quite satisfactory. The re- 
maining difference has to be ascribed to contributions of 
higher order in A. Finally, inserting Eqs. ([64]) and (|65|) 



Delta=0.3 
Delta=0.5 ■ 
Delta=0.9 ■ 




0.5 



1.5 
T 



2.5 



to first order in A, therefore from Eq. (|27|) 

^)(T) = A 2 +4A^ lc$ +1 (T). 



1=1 



Inserting the amplitudes from Eq. (I61| . we find 



FIG. 



(Color online) Functions S (T)/A 2 and 



Im(5 e (T)/A 2 ) for the linear rise model with site 1 ini- 
(63) tially (t = 0) occupied, and the rest empty. Various dashed 
lines: exact numerical calculations for A = 0.3, 0.5 and 0.9 
according to key. Adjacent continuous lines: expressions of 
Eqs. and 



A 2 



S«(T) = A 2 +4^^H^ 



i=i 



y2(i+s) + l 



s=0 



{2(l + s) + l)(l + s)s\(2l + s)\ 



A 2 

— sin 2T, 
2T 



(64) 



where the steps to perform the double sum are detailed 
in the Appendix. 



into Eq. (|26[) we find that to lowest order in A for the 
S (T) and S e (T), but keeping A to all orders in the later 
time-evolution, 



(*) = 



sinT 



T 



0f(i_T)3/2(l +a 2)2 
7T 
4 



sin(2(i-T) 



T 



(66) 



The beauty of this expression is that it shows the effect 
of a finite rise-time T, very cleanly: 1) it adds a phase 
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T to the argument of the sine, and 2) it adds a factor 
(sin T)/T to the amplitude of the sudden approximation 
(the T — > limit.) One sees that for linear switching a 
finite T always reduces the amplitude. Repeating the 
same steps for the c„, starting from Eq. (|20|) . we arrive 
at 

c W m _ »-> + a 2 (n-2)]A 

" W ^ ~ T)3/ 2 (l + « 2 ) 2 

^Bin^-T)-(n-i)|+r), (67) 

which shows the same memory effects. Previously, the 
reduction or enhancement of the amplitude had been dis- 
cussed only for adiabatic switching: see Robinson, [§] and 
references therein. Robinson studied the model of a dis- 
crete state coupled to a continuum extending from a finite 
threshold to infinity. The explicit forms of the continuum 
state density and of the adiabatic switching-on function 
for discrete to continuum coupling, f(t), were left un- 
specified. Making only very general assumptions he was 
able to show that the asymptotic decay was controlled 
by the Fourier transform of the time-derivative f(t), but 
no results for specific attenuation functions were given. 
In this paper the continuum is that of the tight binding 
model; it has finite width, and Robinson's methods do 
not apply. By working with the Fano-Anderson model, 
we have obtained explicit forms, Eqs. (|66|) and ([67]), for 
the memory effects. 



law and that the period of the oscillations, according to 
Eq. ((SB]) or (gSJ), is tt. Neither result is affected by the 
finite rise time T; only the envelope and phase of the 
oscillations are affected by the initial-switching process. 
Finally, we have also plotted the prediction from the ap- 
proximate expression, eq. 1661 On the scale of the figure 
it cannot be distinguished from the green line previously 
obtained, thus confirming the accuracy of our approxi- 
mate analytic expressions. Similarly, we have tested the 




170 172 174 176 178 180 

time(t) 



FIG. 8: (Color online) Survival probability in the Linear rise 
model, for T — 1 and A = 0.3. The continuous (red) line is 
the exact result. Practically hidden by it, the (green) dashed 
line is the value from the asymptotic expressions in the text. 
For comparison, the sudden approximation is also shown, as 
a (blue) dotted line. 



IV. RESULTS FOR LINEAR SWITCHING 

Amplitudes at t — T. To illustrate our results we 
choose A = 0.3 as in Fig.l. As has been shown, Eq. 
(l32f . the lifetime is r = I/7 = 5.3. We only consider 
values of the rise time T«t; otherwise the rise and de- 
cay regimes would be mixed. We first computed the oc- 
cupation amplitudes, c„(T), by numerical solution of the 
time dependent Schrodinger equation as described in [27l ] . 
These were compared to the approximate analytical ex- 
pressions, eq. (|6ip . derived in Section III. The agreement 
with the exact values of c n {T) is excellent, moreso for 
small values of T. But even at T = 2 the absolute values 
are very close: in the format "exact (approximation)" , 
we have for n = 2: 0.222(0.219), n = 3: 0.136(0.136), 
n = 4: 0.0676 (0.0675), n = 5: 0.0275 (0.0274). 

Asymptotic survival probabilities. To compute the 
asymptotic amplitude (t). we then use Eqs. (f2"o) . 
(l27f and the above exact values for the amplitudes at 
t = T . This is compared to the exact numerical solution 
of the TDSE in Fig. El for the case of rise time T = 1. 
Qualitatively similar oscillations are found for other T, 
even as high as T = 2. In all cases the prediction fom Eq. 
(|28p is excellent. Note also that there is a clear difference 
from the prediction given by the "sudden approximation" 
(rise time T — 0.) We stress that the envelope of the os- 
cillations in \ci (t)\ 2 follows the predicted t~ 3 algebraic 



accuracy of Eq. (|6T[) by comparison to exact TDSE cal- 
culations: The agreement found for the other cli\t) is 
similar to that for the amplitudes of the first site. Fur- 
ther, we have found that the additional phase T in the 
argument of the sine, in Eq. f|66[) has the same value 
whatever the shape of the switching function, providing 
that the average value of A v (t) is half of the ultimate 
value A, over the interval < t < T. Finally, we have ex- 
plored the accuracy of the ansatz c n (t) = effl (t) + effl (t) 
for n = 1,2,...: For n = 1 and the range of times of 
Fig. 1, the exact (numerical) result and the values from 
the ansatz are indistinguishable on the scale of that fig- 
ure. Similar agreement is found for other small values of 
n. Support for the validity of this ansatz, for the sud- 
den approximation, can be found in the approach based 
on Laplace transforms presented in [301 ] : the asymptotic 
term is an approximation to the contribution due to the 
Bromwich paths, and to this should be added that of the 
poles, which produce the exponential decay. 

Exponential decay regime In Fig. [9] we compare the 
exact survival probability |ci| 2 , again determined from a 
numerical solution of the TDSE, to the predictions from 
Eq. (|40[) . Three representative values of T are shown. 
Except at small times t < T, the exponential contri- 
bution dominates and accurately follows the behaviour 
of the survival probability. Even for values of T well 
beyond those suggested by the "sudden approximation" 
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condition (TA < 1) [29| one sees that the decay is expo- 
nential with the same decay constant. 




o I > . . . . > — 

2 4 6 8 10 12 14 

time ( t ) 



FIG. 9: (Color online) Linear rise model, A = 0.3: Survival 
probabilities for T = 0.5, 1 and 2 from left to right. Contin- 
uous (red) lines: exact. Dotted (blue) lines: the exponential 
decay component, Eq. (|40p . 




5 10 15 20 25 30 35 40 45 
time ( t ) 



FIG. 10: (Color online) |ci(t)| 2 for, from top to bottom, rise 
times T = 0.5, 1 and 2. Linear rise model with A = 0.3, 
but now site 2 is initially occupied, and the others empty. 
Continuous (red) lines: exact. Dotted (blue) straight lines: 
the exponential decay component, Eq. (I40[) . 

In Figs. [5] and [5] we assumed that at t — only the first 
site was occupied. But, as already stressed, Eqs. (l4Tfl) are 
valid for any set of occupation amplitudes at t = 0. As a 
second example, we explore what occurs when only site 
2 is occupied at t = 0. Fig. [TU] shows the occupation 
probability |ci| 2 for the linear rise model and several val- 
ues of T. Exact (numerical) solutions of the TDSE are 
compared to values from Eq. (|40| . (The relevant set of 
"d q ,q = 1,2,.. was generated by solving the TDSE us- 
ing Ai(t).) The exponentially decaying contribution is 
dominant, except at very small times during which the 
amplitude c±(t) builds up. The maximum occupation is 
of order A 2 , and is lower for larger T, because site 2 is de- 
caying outwards during the switching-on process before 
site 1 can begin to fill. Superimposed on the exponential 



there are pronounced oscillations, due to terms neglected 
when only the pole at z g — ~ia was retained. The aver- 
age trend is well explained by the exponential. 

We have made similar comparisons when initially only 
site n = 5, 10 or 20 is occupied. In all cases the trend is 
similar: there is an initial period where the occupation of 
site 1 accumulates and its subsequent exponential decay 
is well described by Eq. (jlO"]) . Increasing the value of A 
increases the decay constant, so the changeover from ex- 
ponential to the asymptotic decay regime occurs sooner. 



V. SUMMARY AND DISCUSSION 

In the context of Longhi's version of the Fano- 
Anderson model we have obtained explicit expressions 
for so-called "memory" effects. In other models these 
can be studied only numerically or discussed in very 
general terms, leading to few detailed predictions. The 
switching-on process occurs during a time interval T be- 
tween the original preparation of the system in site 1, 
and the stabilization of the coupling to the external sites 
at its final value gA. We have shown that no matter 
what profile A„ (t) is chosen, the site amplitudes oscillate 
asymptotically with period 7r, with an envelope having 
its leading term proportional to i~ 3 / 2 . Furthermore, we 
have also shown that when exponential decay occurs, the 
decay constant is the same for all site amplitudes. 

For the special case of linear switching, we have de- 
rived series expansions for the contributions linear in A 
to the site amplitudes at time T, Eq. (|61[) . We have 
found that to this order the asymptotic amplitude of site 
1 is given by Eq. (|66[) . which shows in a very clean way 
that memory effects introduce a delay T in the phase, 
and that the modulation of the amplitude is given by a 
factor (sin T)/T. Combined with the exponentially de- 
caying contribution they allow accurate predictions in the 
intermediate range of times where interference between 
the two contributions occurs, and this opens the way to 
examine memory effects in this range of times. 

The methods of Section III, could be applied to 
quadratic or other forms of initial switching, allowing a 
systematic study of memory effects. Our predictions for 
the occupation amplitudes in the initial stages should be 
amenable to experimental verification using the waveg- 
uide analogue of the tight-binding model. 
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APPENDIX A 



To sum the double series in Eq. (|64|) we make a change 
of index to k — I + s 



rp2(l + S ) + l 



(2(Z + s) + l)(Z + s)s!(2Z + s)! 
(*-*) 



ySfe+l fc ; ;. . ,2 



E 



( fc( 2 fc+l) ^ s! ( 2 fc- s )! ! 
and to perform the sum over s, we write 

' (k-s) 2 ^ (fc-s) 2 ^ (k-s) 2 



(Al) 



E 

s=0 



E 



s\{2k-s)\ ^s!(2fc-s)! < s\(2k-s)\ 



E 



E 



(fc- S ) 2 



E 



(jfe - 2k + s')' 2 



^sl(2k-s)\ ^ o (2k-s')\(s')\ 

A (fc-s) 2 _ ^ (fc-s) 2 
^s!(2fc-s)! ^s!(2fc-s)! 1 ' J 

s— v ' s— v 7 



where we have introduced s' = 2k — s to rewrite the 
second sum. Next, noticing that 



The double series in Eq. (|65|) is summed using the 
same methods: We again write k = I + s to replace the 
sums over I and s by sums over k and s: 



Qd) = A 2 y H fc ^ y(2fc-2 8 -l) 2 



^fc(2fc-l)^ S !(2fc- S -l)! 



(A7) 



and to evaluate the sum over s we proceed as before 



(2fc - 1 - 2s) 2 
Z. S !(2A;-1- S )! 

2fe-l 

s!(2fc-l-s)! ~ ^ s!(2fc - 1 - s)\ 
-'' 1 (2k - 1 - 2s) 2 {2k - 1 - 2(2fe - 1 - s')) 



(2fc - 1 - 2s) 2 ^ (2fc - 1 - 2s) 2 



E 



- s!(2fc-l-s) 



E 



'112 



! ( S ')!( 2 fc-l-s')! 



' A 1 (2fc-l-2s) 2 ^ (2fc-l-2s) 2 



E 



^ a l(2fc-l-*)! ^s!(2fc-l-s)! 



E 



(A8) 



and therefore 



' {k- S f ^ (fc_ s )2 



E 



^ s!(2fc _ s)! ^ s!(2 fc - s)! 



E 



we finally arrive at 

k 



V (fc-^) 2 (fc-s) 2 

Z^ c !C9i. - ell 9^ 



_ o s!(2fc-s)! 2^s!(2fc- S )! 



(A3) 



(A4) 



fc-i 

E 



(2k - 1 - 2s) 2 _ 1 ^ (2k - 1 - 2s) 2 



_. s!(2 fc-l-s)! 2^s!(2fc-l-s) 



E 



s=0 



(A9) 



Next we decompose the numerator as 
(2k - 1 - 2s) 2 = (2k - l) 2 - 8(k - l)s + 4s(s - 1) (A10) 



Using the identities 



2A- 



s=0 
2k 



2/,- 



]Tsf 2fc ) = 2 2 R 



^s(s-l)( 2fc ) = 2 a *- 1 *(2fc-l) (A5) 



s=0 



and decomposing the numerator of Eq. (|A3[) as (fc— s) 2 = 
fc 2 - (2k - l)s + s(s - 1), Eq. ([64]) simplifies to 

^(T ) = ^EH fc K = A 2 ^.(A6) 



2T ^ 7 (2fc + 1)! 



2T 



and use Eq. (|A5|) . with replacement 2k — > 2fc — 1 to 
rewrite the sum as 



(2fc-l-2s) 2 2 2fe 



2fc-2 



fe-1 

E s i(2fc-l- s )! (2fc-2)f 



We now insert this in Eq. (|A7j) and find 



A 2 g ( _ )fc (2T) 2fe 



2T 



k=l 



(2k)\ 



(All) 



^(l-cos2T) = ^sin 2 T. (A12) 
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